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图 象 作 为 更 直接 更 丰富 的 信息 载体 ， 正 在 成 为 越 来 越 重 要 的 研究 对 象 。 昆 虫 是 生物 类 群 中 物种 最 
丰富 、 数 量 最 多 的 一 类 ， 因 此 也 是 最 为 重要 的 生物 类 群 之 一 。 研 究 昆虫 图 象 的 处 理 和 分 析 对 昆虫 形态 
学 、 昆 虫 生态 学 以 及 昆虫 测报 所 需 的 田间 抽样 调查 具有 重要 意义 ， 最 终 将 实现 田间 昆虫 数量 调查 自动 
计数 的 计算 机 视觉 系统 。 对 昆虫 图 像 的 处 理 和 识别 来 说 ， 研 究 的 兴趣 在 于 图 象 中 的 昆虫 部 分 ， 即 昆 时 
部 分 作为 前 景 。 为 了 识别 和 分 析 目 标 ， 需 要 将 这 些 有 关 区 域 分 离 鉴别 出 来 ， 即 将 图 象 中 昆虫 部 分 与 背 
景 分 割 开 来 。 而 图 象 分 割 的 重要 步骤 之 一 就 是 边缘 检测 。 


1 材料 与 方法 
1.1 材料 

本 文采 用 在 世界 范围 内 广泛 分 布 的 麦 长 管 蚜 Macrosiphum avenae (Fabricius) 无 翅 成 蚜 的 灰 度 图 象 
(356 X 591pixel， 为 便于 排版 ， 本 文中 图 象 均 按 比例 进行 了 缩小 ) 进行 边缘 检测 的 研究 。 边 缘 检 测算 
法 在 Borland C++S. 0 环境 实现 ， 软 件 运 行 环境 为 [BM PC 兼容 机 (Intel Pentium 133M，16M 内 存 )。 


采用 的 边缘 检测 算法 包括 : 梯度 算 子 ，Laplacian 算 子 ，Sobel 算 子 ，Kirsch 算 子 ，Marr-Hildrith 算法 ， 
Canny 算法 ，Shen-Castan 算法 ， 以 及 数学 形态 学 算 子 。 


1.2 边缘 检测 概述 

在 图 象 处 理 和 分 析 中 ， 边 缘 检 测 (Edge Detection) 的 应 用 最 为 普遍， 描述 边缘 检测 或 边缘 增强 
(Edge Enhancement) 的 算法 最 多 41。 昆虫 图 象 中 边缘 形成 昆虫 的 四 周 轮廓 ， 是 昆虫 与 背景 的 边界 。 
如 果 昆 虫 图 象 中 虫 体 的 边缘 可 以 精确 的 鉴别 ， 则 整个 昆虫 虫 体 所 在 的 区 域 就 可 以 确定 下 来 ， 从 而 虫 体 
的 一 些 基本 图 象 特征 如 形状 、 面 积 、 周 长 等 可 以 得 到 测量 。 边 缘 检 测 是 图 象 分 割 过 程 的 一 部 分 5:31， 
一 旦 确定 昆虫 图 象 中 昆虫 的 边界 ， 就 可 将 昆虫 虫 体 与 背景 分 割 ， 得 到 昆虫 图 复 ， 进 而 可 对 昆虫 图 象 进 
行 特征 提取 等 操作 。 


1.3” 几 种 边缘 检测 方法 简介 
在 灰 度 图 象 中 ， 边 缘 就 是 在 特定 位 置 出 现 的 灰 度 变化 ， 如 昆虫 图 象 中 虫 体 与 背景 之 间 具 有 不 同 的 

灰 度 ， 从 而 形成 边缘 ， 即 边缘 是 灰 度 不 连续 的 结果 。 灰 度 变 化 越 大 ， 边 缘 越 容易 检测 !L3]。 常 用 的 边 

缘 检测 算法 有 三 类 ， 即 微分 算 子 、 模 板 匹 配 法 和 基于 边缘 模型 的 算法 [1。 

1.3.1 微分 算 子 : 若 边 缘 定义 为 灰 度 的 阶 跃 变 化 ， 对 灰 度 变化 敏感 的 运算 就 可 以 作为 边缘 检测 算 子 。 

导数 实质 上 反映 了 变化 率 的 大 小 ， 因 此 可 以 对 图 象 求 导 检测 边缘 。 最 简单 的 导数 算 子 是 梯度 算 子 ， 它 

反映 了 边缘 的 强度 信息 (或 称 为 边缘 啊 应 ，Edge Response) 和 边缘 象 素 的 方向 信息 。 
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基于 二 阶 导数 的 检测 方法 称 为 拉 普 拉 斯 (Laplacian) 算 子 。 在 灰 度 边缘 很 了 尘 时 ， 梯 度 算 子 较 有 
效 ， 但 对 于 边缘 过 渡 区 较 宽 时 ， 用 二 阶 偏 导数 更 有 优势 。 在 检测 结果 中 ， 边 缘 处 于 最 大 正 值 和 最 小 负 
值 之 间 的 过 零点 (Zero Crossing)。 
1.3.2 模板 匹配 法 : 模板 匹配 法 是 应 用 离散 的 小 模板 作为 边缘 的 模型 ， 而 不 是 直接 应 用 微分 算 子 。 
模板 匹配 法 既 可 以 是 用 模板 描述 边缘 ， 对 边缘 进行 检测 ， 也 可 以 是 对 微分 算 子 的 一 种 近似 。 而 后 者 是 
更 为 常用 的 方法 上 ,31。 这 里 采用 了 两 种 最 为 常用 的 模板 匹配 算法 : Sobel 算法 和 Kirsch 算法 [1， 因 为 它 
们 能 给 出 较 好 的 检测 结果 。 
1.3.3 ”基于 边缘 模型 的 算法 : Marr Hildreth 算法 。Marr (1980) 试图 将 生物 视觉 的 理论 应 用 于 机 器 
视觉 ， 提 出 了 Marr-Hildreth 算 子 1。 该 算 子 在 使 用 Laplacian 算 子 的 则 时 ， 考 虑 了 边缘 宽度 和 图 象 噪 
声 等 问题 ， 把 Gaussian 滤波 和 Laplacian 滤波 器 结合 起 来 ， 产 生 不 同 大 小 的 卷 积 模板 从 而 检测 不 同 尺 度 
的 边缘 。 具 有 宽 的 Gaussian 滤波 器 检测 宽 的 、 绥 变 的 边缘 ， 罕 的 Gaussian 滤波 器 检测 窗 的 、 尖 锐 的 各 
细节 的 边缘 。(2) Canny 边缘 检测 器 。Canny (1986) 定义 了 一 组 边缘 检测 器 应 符合 的 条 件 ， 并 描述 
了 实现 这 些 目标 的 最 优 方 法 ， 称 为 Canny 方法 561。Canny 认为 边缘 检测 器 应 是 一 个 平滑 噪声 、 定 位 边 
界 的 卷 积 滤波 器 。Canny 对 该 方法 提出 的 条 件 可 以 描述 为 : 由 算法 应 只 对 边缘 响应 ， 且 应 能 检测 到 所 
有 的 边缘 ， 即 错误 率 最 小 。 四 检测 到 的 边缘 象 素 与 实际 象 素 间 的 距离 应 尽 可 能 的 小 。 加 在 只 存在 一 个 
边缘 的 部 位 不 应 检测 到 多 个 边缘 。 (3) Shen-Castan 方法 。Shen 及 Castan (1992) 对 边缘 检测 器 的 一 
般 形 式 看 法 与 Canny 相同 ， 即 先 对 噪声 进行 平滑 ， 再 进行 卷 积 ， 然 后 对 边缘 象 素 进 行 搜索 、 定 位 口 。 
该 滤波 器 使 用 的 平滑 滤波 器 是 无 了 眼 对 称 指数 滤波 器 (Infinite Symmetric Exponential Filter，ISEF )。 
Shen-Castan 还 采用 了 假 过 零点 消除 (False Zero Crossing Suppression) 方法 ， 使 边缘 的 定位 更 准确 。 
1.3.4 数学 形态 学 算法 : 数学 形态 学 是 基本 运算 为 腐蚀 (erosion) 和 膨胀 《dilation)， 由 此 可 以 衍生 
出 其 他 运算 如 开 (opening)、 闭 (closing)、 击 中 (hitting)、 击 不 中 (missing) 等 。 这 些 运算 都 是 运 
用 某 种 称 为 结构 元 素 (structuring element) 的 模板 进行 类 似 卷 积 的 运算 实现 的 ， 实 际 上 ， 结 构 元 素 的 
选择 往往 决定 了 形态 学 运算 的 性 质 '8 "1 。 数 学 形态 学 的 边缘 检测 实际 上 是 对 原 图 象 进 行 腐蚀 运算 ， 得 
到 一 个 新 的 图 象 ， 然 后 求 原 图 象 与 新 图 像 的 差 ， 得 到 边缘 图 象 。 本 文中 形态 学 边缘 检测 是 对 二 值 图 和 象 
进行 的 。 


2 ”结果 与 分 析 


应 用 上 述 边缘 检测 算法 对 麦 长 管 蚜 图 象 (图 1: a) 进行 检测 ， 结 果 见 (图 1: b 一 i)。 

梯度 算 子 由 于 只 是 对 图 象 一 阶 导 数 的 差分 近似 ， 只 能 对 灰 度 变化 剧烈 、 即 尖锐 的 边缘 有 了 响应， 而 
对 于 灰 度 变化 缓慢 ， 平 缓 的 边缘 响应 很 弱 或 没有 响应 ， 因 此 ， 应 用 梯度 算 子 对 昆虫 图 象 的 检测 结果 
(图 1: b) 显示 为 离散 的 点 形成 的 模糊 的 蚜虫 体 轮廓 。 但 如 果 没 有 原 图 (图 1: a) 作 参 照 ， 您 怕 难 以 
鉴别 出 该 图 象 为 蚜虫 的 边缘 图 象 。 

拉 普 拉 斯 算 子 是 图 象 的 二 阶 导 数 ， 由 于 二 阶 导 数 的 性 质 ， 使 拉 普 拉 斯 算 子 对 缓 变 的 、 宽 的 边缘 也 
能 有 较 强 的 响应 ， 因 此 拉 普 拉 斯 算 子 的 检测 结果 (图 1: c) 显示 了 较为 连续 的 虫 体 轮 廓 ， 触 角 、 腹 管 
也 有 了 轮廓 。 但 用 于 进一步 的 图 象 分 割 和 特征 提取 还 达 不 到 要 求 。 

在 灰 度 过 渡 比 较 尖 锐 且 曲 声 较 小 时 ，Sobel 算 子 效果 较 好 。 但 Sobel 算 子 仅 是 对 梯度 算 子 的 改进 ， 
其 检测 结果 (图 1: d) 与 梯度 算 子 相 比 ， 略 为 改善 ， 对 一 些 细 节 的 检测 上 (如 触角 、 腹 管 )，Sobel 算 
子 不 如 拉 普 拉 斯 算 子 。 图 1: e 显 示 了 Kirsch 算 子 的 检测 结果 。Kirsch 算 子 虽然 采用 8 个 模板 对 图 象 进 
行 卷 积 ， 但 对 蚜虫 图 象 的 检测 结果 与 Sobel 算 子 差别 不 大 。 

Marr-Hildreth 算法 对 边缘 非常 敏感 。 通 过 分 析 比 较 ， 对 于 蚜虫 图 象 (图 1: a) 的 边缘 检测 ， 我 们 
取 o=3.0， 采 用 15x15 和 27x27 的 模板 进行 两 次 高 斯 平滑 ， 并 分 别 进 行 LOG 卷 积 ， 求 取 过 零点 作 
为 边缘 象 素 ， 将 两 次 计算 得 到 的 边缘 象 素 进行 合并 ， 与 上 述 算法 比较 ， 得 到 较 好 的 边缘 图 象 (图 1: 
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f)。 但 由 于 模板 较 大 ， 使 得 计算 过 程 较 费 时 ， 而 且 该 算法 在 不 含 边缘 的 区 域 检测 到 了 边缘 ， 这 是 我 们 
所 不 希望 的 。Marr-Hildreth 算法 能 够 检测 到 较 完整 的 边缘 ， 但 也 检测 到 了 不 存在 的 边缘 ， 使 得 图 象 分 
割 增加 了 难度 。 对 Canny 算法 ， 经 过 比较 分 析 ， 标 准 差 参数 取 值 为 s= 0.4 可 取得 理想 的 边缘 图 象 。 
由 图 1. g 可 看 出 ，Canny 方法 的 检测 结果 较 Marr-Hildreth 方法 理想 ， 没 有 包含 不 应 有 的 边缘 ， 而 且 
虫 体 边缘 较为 连续 ， 许 多 细微 部 位 如 腹 管 边缘 清晰 可 见 。 图 1. h 说 明 Shen-Castan 方法 能 够 检测 出 非 
常 细微 的 边缘 ， 轮 廊 较 完整 ， 虫 体 的 整体 边缘 清晰 而 简单 。 但 也 检测 到 了 一 些 不 存在 的 边缘 〈 图 象 左 
上 部 )， 这 可 能 是 由 于 Shen-Castan 方法 没有 采用 Canny 方法 中 的 多 啊 应 限制 (Multiple-Response Con- 
straint) 指标 ， 因 此 它 的 算法 可 能 会 对 噪声 或 模糊 边缘 产生 假 的 边缘 响应 。Canny 算法 和 Shen-Castan 
算法 检测 到 的 边缘 都 能 用 来 进行 进一步 的 图 象 处 理 ， 如 图 像 分 割 。 





图 1 以 不 同 边 缘 检 测算 法 对 麦 长 管 蚜 Macrosiphum avenae (Fabricius) 图 象 的 处 理 结果 
a. 麦 长 管 蚜 无 翅 成 蚜 图 象 ; b. 梯度 算 子 检测 到 的 边缘 图 象 ; c. 拉 普 拉 斯 算 子 (Laplacian) 检测 到 的 边缘 图 象 ; 
d. Sobel 算 子 的 检测 结果 ; e. Kirsch 算法 的 检测 结果 ; f. Marr - Hildreth 算法 的 检测 结果 ; g. Canny 算法 的 
检测 结果 ; h. Shen - Castan 算法 的 检测 结果 ; i. 形态 学 边缘 提取 算法 的 结果 
Fig. 1 The detection result images of Macrosiphum avenae (Fabricius) by using various edge detectors. 
a. Image of M. avenae; b. Detection result of gradient operator; c. Detection result of Laplacian operator; 
d. Detection result of Sobel operator; e. Detection result of Kirsch operator; f detection result of 
Marr-Hildreth algorithm; g. Detection result of Canny algorithm; h. Detection result of 


Shen-Castan algorithm; 1. Detection result of morphological method 


在 蚜虫 图 象 的 边缘 检测 中 ， 数 学 形态 学 算法 的 确 表现 了 其 先进 性 和 准确 性 。 由 图 1. i 可 以 看 出 ， 
检测 的 结果 中 蚜虫 图 象 连续 且 清 晰 ， 而 且 非 常 简 炼 。 而 该 算法 的 运算 速度 也 是 最 快 的 。 


3 讨论 


本 文 尝试 应 用 不 同 的 边缘 检测 算法 进行 昆虫 图 象 的 边缘 检测 研究 ， 结 果 表 明 ， 微 分 算 子 和 模板 匹 
配 算法 由 于 检测 结果 模糊 且 不 连续 ， 均 无 法 满足 蚜虫 图 象 分 割 的 要 求 ; Marr-Hildreth 方法 对 边缘 很 敏 


208 昆 虫 学 报 41 卷 


感 ， 能 够 检测 出 很 显著 的 边缘 ， 但 也 检测 到 不 存在 的 假 边 缘 。 比 较 理 想 的 算法 是 Canny 方法 、Shen- 
Castan 方法 和 数学 形态 学 边缘 检测 算法 ， 这 三 种 算法 检测 到 的 边缘 都 可 进一步 用 于 图 象 分 割 。 但 相 比 
较 而 言 ， 数 学 形态 学 算法 因 其 算法 的 先进 性 、 检 测 结果 的 近乎 完美 和 运算 速度 快 等 优势 ， 将 作为 图 象 
分 割 过 程 中 的 首选 边缘 检测 算法 。 而 Canny 和 Shen-Castan 方法 算法 相 比 则 Canny 方法 对 蚜虫 图 象 边 
缘 更 为 敏感 ， 其 检测 结果 亦 可 用 于 图 象 分 割 。 基 于 模型 的 三 种 算法 都 可 以 用 不 同 的 参数 值 来 调整 检测 
右 的 检测 行为 ， 但 Shen-Castan 方法 似乎 对 参数 变化 的 响应 不 明显 ， 在 不 同 参 数值 下 检测 结果 没有 区 
别 ; 而 Marr-Hildreth 算法 和 Canny 算法 对 参数 值 的 变化 反应 明显 ， 本 文 所 选 的 检测 结果 均 取 其 不 同 参 
数值 下 的 最 优 检测 结果 图 象 。 

由 于 昆虫 体形 变异 很 大 ， 体 色 复 杂 多 变 ， 既 使 是 昆虫 的 灰 度 图 象 也 很 复杂 ， 因 此 昆虫 图 象 的 边缘 
检测 就 显得 很 困难 ， 对 检测 算法 的 选择 和 验证 也 就 很 重要 了 。 对 蚜虫 而 言 ， 由 于 虫 体 非 常 小 ， 再 加 上 
昆虫 是 活动 的 ， 这 不 仅 增 加 了 蚜虫 图 象 采样 的 困难 ， 而 且 量 化 的 图 象 往往 也 很 难 达 到 令 人 满意 的 程 
度 ， 而 图 象 的 质量 对 边缘 检测 来 说 也 是 很 重要 的 ， 因 此 蚜虫 图 象 的 处 理 和 分 析 尤 为 困难 。 另 外 ， 蚜 贝 
与 其 寄主 植物 的 密切 关系 也 增加 了 蚜虫 图 象 处 理 的 难度 。 幼 小 的 蚜虫 总 是 生活 在 寄主 植物 叶片 上 的 ， 
因此 多 数 蚜 虫 图 象 都 是 以 植物 叶片 为 背景 ， 而 长 期 的 自然 选择 和 进化 使 蚜虫 具有 保护 色 ， 虫 体 大 多 数 
时 间 呈 黄 绿 色 ， 与 植物 叶片 很 难 区 分 ， 反 映 在 灰 度 图 象 上 则 是 灰 度 值 的 变化 不 大 ， 边 缘 过 渡 缓 慢 。 

除了 对 边缘 检测 算法 进行 选择 和 验证 外 ， 对 所 选择 的 算法 结合 昆虫 特点 进行 改进 也 是 必要 的 。 另 
外 ， 检 测 到 的 边缘 不 完全 连续 ， 还 无 法 直接 用 来 对 虫 体 和 背景 进行 分 割 ， 还 需要 进行 边界 闭合 运算 ， 
使 检测 到 的 边缘 成 为 完全 连续 的 边界 ， 这 样 进行 分 割 就 容易 进行 了 。 
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